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STRUCTURAL PARAMETERS FOR GLOBULAR CLUSTERS IN M31 

SongWang,'-2'3 JunMa'-'' 
A], in press 

ABSTRACT 

In this paper, we present surface brightness profiles for 79 globular clusters in M31, using images observed 
with Hubble Space Telescope, some of which are from new observations. The structural and dynamical param- 
eters are derived from fitting the profiles to several different models for the first time. The results show that in 
the majority of cases. King models fit the M3 1 clusters as well as Wilson models, and better than Sersic models. 
However, there are 1 1 clusters best fitted by Sersic models with the Sersic index n > 2, meaning that they have 
cuspy central density profiles. These clusters may be the well-known core-collapsed candidates. There is a 
bimodality in the size distribution of M3 1 clusters at large radii, which is different from their Galactic coun- 
terparts. In general, the properties of clusters in M31 and the Milky Way fall in the same regions of parameter 
spaces. The tight correlations of cluster properties indicate a "fundamental plane" for clusters, which reflects 
some universal physical conditions and processes operating at the epoch of cluster formation. 
Subject headings: galaxies: individual (M31) - galaxies: star clusters general - galaxies: stellar content 



1. INTRODUCTION 

It is well known that, studying the spatial structures and 
dynamics of globular clusters (GCs) is of great importance 
for understanding both their formation condition and dy- 
namical evolution with in the environment of their galaxies 
dMclaughlin et al.l2008l) . For example, these clusters are ideal 
laboratories for detailed studies on two-body relaxation, mass 
segregation, stellar collisions and mergers, and core collapse 
(^eylan & Heggie 1997). The correlations of structures with 
galactocentric distance can provide information on the role of 
the galaxy tides towards the clusters, while the distribution 
of ellipticity can shed light on the primary factor for cluster 
elongation. In addition, comparisons of structures of GCs lo- 
cated in different environment of galaxies offer clues to differ- 
ences in the early formation and evolution of the galaxies or 
in their subsequent accretion histories (Bellazzini et al. 2003; 
iMackev et alj|2007h . The "fundamental plane" for clusters 
in parameter space reflects universal cluster f ormation con- 
diti ons, regardless of their host environments (iBarmby et alJ 
l2009h . 

The structural and dynamical parameters of clusters are of- 
ten determined by fitting the surface brightness profiles to 
structure models, combined with mass-to-light ratios esti- 
mated from velocity dispersions or population-synthesis mod- 
els. An accurate and well resolved density profile can be 
obtained by studying t he distribution of int egrated light cou- 
pling with star counts (Federici et al. 2007). Several models 
are often used in the fits: the empirical, single-mass, mod- 
ified isothermal spheres (King 1962, 1966; Wilson 19751); 
the isotropic multi-mass models (Da Costa & Freeman 19761); 
t he anisotrop i c mul ti-mass modelslOunn & Griffin 119791 ; 

flevlah 1988', the power-law surface-density profiles 

Jersic 1968; Elson, Fafl & Freeman 1987). 
The nearest large GC system outside the Milky Way (MW) 
is that of M31, with a distance of ~ 784 kpc from us 
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dStanek & Ga r navich|[T998h . It is so close to us that most 
GCs in it can be well resolved w ith Hubble Space Tele- 
scope (HST). ,Battistini et al.l (I1982I) first estimated core radii 
for several clusters in M31, and subsequ e ntly, a number 
of studies (Pritchet & van den Bergh' '1984* 'Crampton et alJ 
1985; Spassova et al. 1988; Bendinelli et al. 1990, 1993F 
Cohen & Freeman 1991; Fusi Pecci et al.ll994tlBarmbv et al] 
2002; Ma et al, 2006.. ,2007, ,20 IZ) focused on the internal 
structures of M3 1 GCs, including the core radius, half-light 
radius, tidal radius, and ellipticity, using the images from large 
ground-based telescopes and HST. Barmby et al. (2007) de- 
rived structural and dynamical parameters for 34 M3 1 GCs, 
and construct a comprehensive catalog of these parameters 
for 93 M3 1 GCs with corre cted versions of those in a previous 
study ( Barmby et al.l2002h . Combined with the structures and 
dynamics for clusters from the MW, Magellanic Clouds, the 
Fornax dwarf spheroidal, and NGC 5128, these authors found 
the GCs have near-universa l structural pro perties, regardless 
of their host environments. iBarmbv et al.l (l2009l ) found that 
bright young clusters in M31 are larger and more concen- 
trated than old ones, and are expected to dissolve within a 
few Gyr and will not survive to become old GCs. With mea- 
surements of structural parameters for 13 extended c lusters 
(ECs) in the halo regions of M31, Hu xor et"aI1 (1201 1) pre- 
sented that the faintest ECs have magnitud es and sizes similar 
to Palomar-type GCs in the MW halo. Wa ng & Mai (120 12h 
measured structures and kinematics for 10 newly discovered 
GCs in the outer halo of M31, and found that they have larger 
ellipticities than most of GCs in M3 1 and the MW, which may 
be due to galaxy tides from satellite galaxies of M3 1 or may 
be related to the merger and accretion history t hat M31 has 
experi e nced. Using the sam e sample clusters in I Wang & Mai 
(l2012h . lTanvir et all (l2012h found that some G Cs in M31 ex- 
hibit cuspy cores which are well described by ISersicI (11968!) 
models. These authors also confirmed the exist of luminous 
and compact globulars at large galactocentric radii of M31, 
with no counterparts found in the MW. The last three stud- 
ies extended the structural analysis of clusters in M3 1 out to 
7?gp ~ 100 kpc, providing important information on the accre- 
tion history of M3 1 outer regions. 
In this paper, we determined structures and kinematics for 
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79 clusters in M3 1 by fitting several structural models to their 
surface brightness profiles. This paper is organized as follows. 
In Section 2, we present the HST observations for the sample 
clusters, and the data-processing steps to derive the surface 
brightness profiles. In Section 3, we determine structures and 
kinematics of the clusters and make some comparisons with 
previous studies. A discussion on the correlations of the de- 
rived parameters is given in Section 4. Finally, we summarize 
our results in Section 5. 

2. DATA AND ANALYSIS PROCEDURES 

2.1. Globular Cluster Sample 

HST has imaged a large fraction of globular clusters (GCs) 
in M31. B armby et al. ( 2002 ) used HST/Wide Field Plan- 
etary Camera 2 (WFPC2) i mages to m easure ellipticities, 
position angles, and best-fit iKingI (1196 6) model (hereafter 
"King model" ) param eters for a large sample of M31 GCs. 
iBarmbv et al.l (120071) determined structures and kinematics 
for M31 GCs by fitting surface brightness profiles to differ- 
ent models, howeve r, for all the GCs but GOOl studied in 
IBarmbv et al.l (120021) . only King model was used. In addi- 
tion, there were some clusters located at the edges of the 
images or observed with only one filter With updated ob- 
servations by HST, new data can be derived for them now. 
So, we decided to re-estimate the structure parameters for 
these GCs in B armby et al. (2002). However, Barmbv et all 
(1^07) determined structur es for GOOl using th ree models 
and new observation, while IBarmbv et al.l (120091) determined 
structure parameters for five clusters (B315, B318, B319, 
B368, and B374) using updated HST data. In addition, there 
is no new observation for B077, which is at the edges of 
the HST images. We also noticed that BH20 and BH21 had 
been classified as stars, while BH24 a galaxy dCaldwell et alj 
l2009h . which would not be discussed again i n this paper. 
The remaining 69 clusters from Barmbv et al. (2002") would 
be included in the sample. Mackev et al. (2007) estimated 
metallicities, distance moduli and reddening values for 10 
newly-discovered halo GCs in M3 1 using the ACS / Wide Field 
Camera (WFC) ini ages. Although IWang & Mai (120121) and 
Tanvir et aP ('2012!) determined structures for t he 10 GCs, 
Wang & Ma (2012) only used King model, while lTanvir et alj 
12012) presented few structure parameters. So, we re- 
estimated the structure parameters for these GCs. Finally, 
there are 79 clusters in our sample. We obtained the com- 
bined drizzled images from the Hubble Legacy Archive. The 
images in the bandpass close to V band (F555W or F606W) 
and / band (F814W) were preferred, otherwise the images of 
F300W, F435W or F475W were selected. The images with 
high resolution were preferentially adopted (WFPC2/PC or 
ACS/WFC), followed by those obtained with WFPC2/WFC. 
Figure 1 shows the spatial distribution for the sample GCs. 

2.2. Surface Brightness Profiles 

The data analysis procedures to measure s urface bright- 
ness profiles of clusters have been described in 'Barmb v et al.l 
(|2b07). When the center positions of these clusters were de- 
termined using the imcentroid task in IRAF, the ellipse task 
was run in two passes to derive the surface brightness pro- 
files. The ELLIPSE showed inability to converge for several in- 
dividual clusters. In thes e cases we first sm oothed the im- 
ages with a boxcar filter (iLarsen et al.l l2002b . and then ran 
the ELLIPSE to derive the density profiles. The overall ellip- 
ticity and position angle (PA) were determined by averaging 
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Fig. 1 . — Location of our sample GCs in relation to M3 1 . The inner el- 
lipse delineates MSl's main stellar disk ((' = 77° and R = 2° ) while the 
oute r ellipse has a radius of 55 kpc and is flattened to hja = 0.6, as given in 
I Richardson et al. (2009). The filled circles and asteiisks represent the sample 
GCs from Barmbv et al. (2002) and Mackev et al. (2007), respectively. The 
two small ellipses near the M3 1 center represent M32 (bottom) and NGC 205 
(top-right). 

the ELLIPSE output in the first pass, with the errors estimated 
as the standard deviation of the mean. Several clusters (B330, 
B468, BH04, BHll, BH29, and NB39) showed eiTors of PA.s 
larger than 15 degrees. We checked the images for these clus- 
ters, and found that the random fluctuations due to individ- 
ual stars (jLarsen et al. 2002) may account for the high errors, 
leading to a more difficult business for accurate measurements 
of PAs Table 1 lists the average ellipticity, PA. and some ad- 
ditional integrated data for the sample clusters. Considering 
that the metallicities for young (< 1 Gyr) and old clusters 
are quite different, first, we averaged the age s for the sam- 
ple 79 clusters from several previous studies (Caldwel l et alj 
12009.. .201 1; Kang et al. 2012 : Wa ng et aL20ia,2012) to dis- 
tinguish young from old clusters. Of which, the re are 13 clus- 
ters (B H23, BH29, NB39, and the 10 GCs froml Mackevetal] 
120071) with no available age values in the literatures, and we 
assumed them as old ones. Finally, there are six young clus- 
ters (B097D, B324, BH05, BH12, DA038, and M091) in our 
sample, all of which are f rom Barmbv et al. (2002). For the 69 
GCs from Barmby et al.' ("2002), the metallicities with uncer- 
tainties from Kang et al. (2012^ were preferentially adopted 
as ou r reference, followed by those of ICaldwell et al.l (120091 
1201 II) for old clusters, while t he solar metal licity was as- 
sumed for young clusters (Barmby et al. 2009). The redden- 
ing values were from Kang et al. ( 2 01 2). while the other in- 
tegrat ed data were from RBC V.5 (iGalleti et al.||2004l l2006l 
120091) . For the 10 GCs from'Mackey et al.' ('2007'), we used 
the integrated data presented in Wang & Ma (2012), includ- 
ing the VI magnitudes, galactocentric distances, reddening 
values, and metallicities (see Wang & Ma 2012, and refer- 
ences therein). In addition, B magnitudes for six of the 10 
GCs are from RBC V.5. Old clusters with no metallicity mea- 
surements are a ssigned with a mea n M31 GC metallicity of 
[Fe/H]= -1.2 ( iHuchraet all 1 1 9911) . while the uncertainties 
of [Fe/H] are assumed to be 0.6 as for the standard devi- 
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ation of the metallic ity distribution of the M31 GC system 
dBarmbv et al.ll2000l) . Clusters with no reddening values are 
assigned with the Galactic reddening in the direction of M3 1 
of E{B - y) = 0.08 ( van den Bergh 1969). 

Raw output from package ellipse is in terms of counts 
s"' pixel"', which needs to multiply by a number (400 for 
ACS/WFC, 100 for WFPC2/WFC, and 494 for WFPC2/PC) 
to convert to counts s"' arcsec"' . A formula was used to trans- 
form counts to surface brightness in magnitude calibrated on 
the VEGAMAG system, 

///mag arcsec"^ - -2.5 log(counts s"' arcsec"') + Zeropoint. 

(1) 

As noted by lBarmby et al.| jlOO T), occasional oversubtrac- 
tion of background during the multidrizzling in the automatic 
reduction pipeline led to "negative" counts in some pixels, 
so we worked in terms of linear intensity instead of surface 
brightness in magnitude. With updated absolute magnitudes 
of the sun Mq (C. Willmer, private communication) listed 
in Table 2, the equation for transforming counts to surface 
brightness in intensity was derived, 

I/Lq pc"^ ^ Conversion Factor x (counts s"' arcsec"'). (2) 

Converting from luminosity density in Lq pc"^ to surface 
brightness in magnitude was done according to 

yLf/mag arcsec"^ - -2.5 log(//L0 pc"^) + Coefficient. (3) 

Table 2 presents the Zeropoints, Conversion Factors, and 
Coefficients used in these transformations for each filter Ta- 
ble 3 gives the final, calibrated intensity profiles for the 79 
clusters but with no extinction corrected. The reported in- 
tensities are calibrated on the vegamag scale. Column 7 
gi ves a flag for each point, which has the same meaning 
as iBarmby etall (1200 7 ) and M claugMin et al.l (120081) defined. 
The points flagged with"OK" are used to constrain the model 
fit, while the points flagged with "DEP" are those that may 
lead to excessive weighting o f the central region s of clusters 
(see lBarmbvlTa l. 2007; Mcla ughUn et al.ll2008l for details). 
In addition, points marked with "BAD" are those individ- 
ual isophotes that deviated strongly from their neighbors or 
showed irregular features, which were deleted by hand. 

2.3. Point-spread Function 

The point-spread function (PSF) models are critical to accu- 
rately meas ure the shapes o f objects in imag es taken with HST 
dRhodes et al. 2006). Bar mbv etal.1 (120021) found that by fit- 
ting models without PSF convolution, the scale radii were sys- 
tematically larger, and the concentrations smaller than those 
estimated from the convolved models. Compared to ground- 
based telescopes, the PSF of HS T is ve ry stable, although it is 
also known to vary with time ( Krist et al. 2011 ). Tiny Tim has 
been the standard modeling software for HST PSF simulation 
for 20 years, with a variety of uses ranging from deconvolu- 
tion, model convolution, PSF fitting photom etry and astrome- 
try, and PSF subtraction (iKrist et al.ll201 Ih . In this paper, we 
derived the ACS/WFC and WFPC2 PSF models using Tiny 
Tim|3, and then the models were fitted using a function of the 
form 

/psF//o = [l+(-R/ro)"]"^^ (4) 

http://tinytim.stsci.edu/cgi-biii/tinytimweb.cgi. 



where ro, a, and /? for each filter are given in Table 4. It 
can be seen that the parameters in Table 4 are slightly dif- 
ferent with those from lBarmbv et al.'('2007') for ACS/WFC in 
the F606W and F814W filters. Barmby et al. (2007) selected 
a number of isolated stars on a number of images, and com- 
bined them to produce a single, average PSF for each filter 
Here we derived a few model PSFs at different positions of 
the camera, and averaged them to produce the final PSF for 
each filter. The discrepancies of these parameters from the 
two studies are due to the different methods, but are negligi- 
ble. The PSF variation over the cluster extent was ignored 
since t he clusters ar e small compared to the camera field of 
view dBarmbv et al.ll2009.) . 

3. MODEL FITTING 
3.1. Structure Models 

We used three structural mode ls to fit star cl uster s urface 
profiles, including King model, ' WilsonI (Il975h . and iSersid 
( 1968 ) model (hereafter "Wilson model" and "Sersic model"). 
Mcla ughlin et al.l (120081) have described the three structural 
models in detail, here we briefly summarized some basic char- 
acteristics for them. 

King model is most commonly used in studies of star clus- 
ters, which is the simple mod el of single-m a ss, is o tropic , 
modified isothermal sphere. iBarmb v etaTI (120071 120091) 
found that M3 1 clusters are better fitted by King models. The 
phase-space distribution function for King model is defined as 

/(£)-{^;^p[-^/'^o]-i^ ^<o, 

where E is the stellar energy, cro is a velocity scale. 

Wilson model is an alternate modified isothermal sphere 
based on the ad hoc stellar distribution function of 
I WilsonI (119751) . These models have more extended enve- 
lope structures than the standard King isothermal spheres 
(Mclaughlin et al. 2008). Several studies presented that Wil- 
son models fit the majority of GCs in the Milky Way (MW) 
and some of its satellit es and NGC 5128 as well as or 
better than King model s dMclaughlin & van der Mareill2005l : 
iMclaughUn et al.l [20081) . The phase-space distribution func- 
tion for Wilson model is defined as 

/(£)oc I ^^Pt-^/'^ol- 1+^/^0' ^<0, (6) 

Sersic model has a R^^" surface-density profile, and has 
been the standard model for parameterizing the surface bright- 
ness profiIes_of_early-tY£egala^^ galax- 
ies (iBaes & Gentild 1201 II) . iTanvir et al.1 (120 12h found that 
some classical GCs in M31 which exhibit cuspy core pro- 
files are well described by Sersic models of index « ~ 2 - 6. 
The clusters with cuspy cores have usu ally b een called post- 
core collapse (see Noyola & Geb hardtl 120061 and references 
therein). The Sersic model is defined as 

m = /o exp[- ln(2) x (R/rof''% (7) 

3.2. Fits 

Before we fitted models to the brightness profiles of the 
sample clusters, the intensity profiles were corrected for ex- 
tinction. Table 2 lists the effective wavelengths of the ACS 
and WFPC2 filters from the Instr ument Handbook. With 
the extinction curve taken from Cardelh et al.l (119891) with 
Rv = 3.1, we derived the values for each filter. 



Wang et al. 




R/arcsec 



Fig. 2. — Surface brightness profiles and model fits to one sample cluster 
B006, with the data of F555W and F814W band from top to bottom. The 
three panels in each line are fits to, from left to right: King model, Wilson 
model, and Sersic model. 

We first convolved the three models with PSFs for the filters 
used. Given a value for the scale radius ro, a dimensionless 
model profile I^od = Imodlh was computed, and then the con- 
volution was carried out. 



mod 



{R\n)) 



Imod(R'/n))hsF [(x - x'), (y - y')] dx'dy', 

where - + y^, and R'^ - y?2 j^^^ appr oxi- 
mated using the equation dU (see lMclaughhn et alj|200 8l for 
details). The best fitting model was derived by calculating 
and minimizing as the sum of squared differences between 
model and observed intensities. 
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VoURi) - toI*^JRi\ro) - /bkg]' 
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in which a background /bkg was also fitted. The uncertainties 
of observed intensities listed in Table 3 were used as weights. 

As an example, we plotted the fitting for one sample clus- 
ter in Figure 2. The observed intensity profile with extinc- 
tion corrected is plotted as a function of logarithmic projected 
radius. The open squares are the data points included in the 
model fitting, while the crosses are points flagged as "P EP" or 
"BAD ", which are not used to constrain the fit ( Wang & Mai 
120121) . The best-fitting models, including the King model, 
Wilson model, and Sersic model are shown with a solid line 
from the left to the right panel, with a fitted /bkg added. The 
dashed lines represent the shapes of the PSFs for the filters 
used. There are some clusters showing individual isophotes 
with ELLIPSE intensities that showed irregul ar features or de- 
viated strongly from their neighbors. As [Mclaughlin et al.l 
(|2008) reported, such bumps and dips may skew the follow- 
ing model fits. In these cases, we first derived the ellipse out- 
put throu gh a boxcar fil ter to make a smoothed cluster profile 
dMclaughhn et al.ll2008h . and then fitted these surface profiles 
using structural models. If some individual isophotes still can- 
not be well fitted, these points were deleted by hand, which 
were masked with "BAD" in Table 3. 

Most profiles of the sample clusters were well fitted by 
the models, except for several clusters with different rea- 
sons. There are one or several bright objects at the inter- 
mediate radii of B097D, GC9, and M091; the shape is very 
loose for DA038; the signal-to-noise ratio (SNR) is low for 
B205; there are one or several bright objects near the outer 
region of B328, B331, and BHll; the images are not well 
resolved for B092, BlOl, B145, BH05, and NB39; sev- 
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Fig. 3. — Images (F555W or F606W) of eight "ring cluster" candidates. 
From the upper left, these are B324, B330, B333, B468, GC3, GC7, BH04, 
and BH29. 

eral clusters (B324, B330, B333, B468, GC3, GC7, BH04, 
and BH29) are lack of a central brightness concentration, 
which may be candidates of "ring clusters". The "ring clus- 
ters" have been reported in the Magellanic Clouds (MCs) 
and M33 (M ackey & Gilmore 2003; Hill & Zaritsky 200^ 
iWerchan & Zaritskvii201 Ic iSan Roman et al.il2012 ), with ir- 
regular profiles such as bumps and dips which may not be 
attributed to the random fluctuations due to a few luminous 
stars. The images of these "ring cluster" candidates in our 
sample are displayed in Figure 3. There are three clusters 
(B018, B114, and B268) showing bumps only in the luminos- 
ity profiles of the F814W, with the intrinsic color of (V-I)q - 
0.95, 1.14, and 1.2, respectively. Some bright redder stars 
may locate at the intermediate radii of these clusters. 

There are 15 clusters (B124, B127, B128, B148, B153, 
B167, B268, B331, BH05, BH12, M091, NB39, B064, B092, 
and B205) with high errors of /bkg- We checked the im- 
ages carefully, and found that several reasons may account for 
those high errors: most of these clusters are not well resolved 
in the images; there are one or several bright objects in the 
intermediate and outer part of clusters B331 and M091; clus- 
ters B064 and B205 only have images in one band (F300W), 
of which the SNRs are low. We should notice that if the fit- 
ted background is too high, the tidal radius may be estimated 
smaller artificially. Only one band of observations can be 
derived for six clusters (B009, B020D, B064, B092, BlOl, 
B205), and none of the bands is close to V band. So, we would 
not include the six objects in the following disc ussions. 

We used the same m et hod m entioned in iBarmbv et al.l 
( I2009h and iWang&Mal ( 120121) to transform the magni- 
tudes from ACS and WFPC2 to V on the vegamag scale 
dHoltzman et al.l l995; Sirianni et al. 2005). For clusters with 
available data of F555W or F606W band, we briefly used 
the extinction-corrected color (V - /)o or (B - V)q to do the 
transformations, while for clusters with no data of F555W or 
F606W band, we first transformed the ACS or WFPC2 mag- 
nitudes to / magnitude using the color (B - /)o, and then com- 
puted the V magnitude using the color {V - I)o- The BVI 
magnitudes and the reddening values are listed in Table 1. 
We estimated a precision of +0.05 mag in the transforma- 
tion, which was propagated through the parameter estimates 
(IBarmbv et aT]|2()07L IWang & Mal l2012[) . 

The mass-to-light ratios {M/L values), which were used to 
derive the dynamical paramet ers, were determined from the 
populati on-synthesis mod els of lBruzual & Chariot (l2003i) , as- 
suming a lChabrie^ (12003b initial mass function. The ages and 
metallicities used to computed M/L values in V band were 
derived as follows. An age of 13 Gyr with an uncertainty of 
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+2 Gyr was adopte d for old clusters, while the derived aver- 
age ages (in Section lZZt were adopted for the young clusters. 
The metallicities with uncertainties are listed in Table 1 . The 
errors for Mj L of old clusters include the uncertainties in age 
and metallicity, while for young clusters, we simply adopted 
10% in MIL as the errors as Barmby et al. (2009) did. 

The basic structures and various derived dynamical param- 
eters of the best-fitting models for each cluster are listed 
in Table 5 to Table 7, with a description of each parame- 
ter/column at the end of each table (see for details of their 
calculation in Mclaughlin et al. 2008). The uncertainties of 
these parameters were estimated by calculating their varia- 
tions in each model that yields x within 1 of the global 
minimum for a cluster, while combined in quadrature with 
the uncertainties in MIL f or the parameters related to it (see 
[Mclau ghlin & v an der Ma rel 2005, for details). 

3.3. Comparison with Previous Determinations 

In this paper, we determined structural parameters for 79 
clusters in M31 by comparing their surface brightness pro- 
files with three structural models. In Figure 4, some estimated 
structural parameters for clusters in this paper were compared 
with t hose from previou s studies (e.g., Barmby et al. 20021 
l2007t IWang&Mal[20T2h . IBarmbv et all (l2002i l2007h pre- 
sented structural and dynami cal parameters in V band for 5 1 
clusters in our sample, while IWang & Mai (120121) determined 
structures for 10 GCs using King model. So we used the re- 
sults on the bandpass close to V band (e.g., F475W, F555W, 
and F606W) and fitted by King model for comparison. It is 
not unexpected to see that most of our par ameters are larger 
than the results from Ba rmby et al.l (|2002[ ^007), since the 
isophotes flagged as "PEP " in this paper may no t be excluded 
from the fitting process in iBarmbv et al.l (120021) . resulting in 
excessive weighting of the inner regions in the fits. In addi- 
tion, the ellipticities presented by Barmby et al. (2002) were 
averaged over diff'erent filters for each cluster, while our re- 
sults in Figure 4 were on the bandpass close to V band. The 
cluster with the largest di screpancy of eUiptic ity is BH05, 
with an estimate of 0. 19 bv lBarmbv et al.1 ( 120021) . and 0.55 in 
this paper As discussed above, the image of BH05 is not well 
resolved, so it is difficult to derive a ccurate ellipticities and 
structural parameters. Ba rmby et al.l (^07) also concluded 
that the shapes of outer parts of faint clusters are strongly af- 
fected by the galaxy background or low SNRs, leading to a 
difficult business to accurately measure the ellipticities. The 
largest scatter of the /?/, is B018, which is ~ 5 pc derived 
by Barniby et al. (2007), while ~ 30 pc in this paper Al- 
though the PSFs and some calibrat ion factors a d opted here 
are slightly different from those of I Wang & Mai (12012b . the 
parameters derived here are in good agreement with their re- 
sults; 

IStrader et al.l (I2009ll201 ll) presented observed velocity dis- 
persions for a number of GCs in M31 using new high- 
resolution spectra from MMT/Hectochelle. These authors 
estimated the M/L values in V band for clusters using the 
virial masses and luminosities, since the virial masses are 
the nominal estimates of the "global" mass of the system 
and are less sen sitive to the accuracy of the King model fit 
(see lStrader et a l. 2009, 2011, for details). Here we estimated 
the MIL values in V band from population-syr ithesis models 
by giving th eir metallicities and various ages (IBarmbv et alj 
120071 120091) . In Figure 5, We presented the compariso n of 
some parameters derived by Barmbv et al. (2007', 12009b and 
this paper with those from .Strader et al., (.201 Ij . There is a 
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Fig. 4. — Comparison of our newly obtained cluster structural parameters 
with previous measurements by Barm bv et al.. (.2002. 2007.) (open circles) 
and . Wang & Mai (20121) (filled circles). 




Fig. 5. — Com paris on of some dynamical parameters from IBarmbv et al] 
120071) (crosses). IBar mbv et al. (2009) (triangles), and this paper (filled cir- 
cles) with those from, Strader et al., (2 01 1.) . 



large discrepancy between the M/L values derived from the 
two methods, and most of the M/L values from population- 
synthesis models are larger than those from observed veloc- 
ity di spersions. The model masses from lBarmbv"etaLl (l2007l 
120091) and this paper are s lightly larger than the virial masses 
from I Strader et al.l (1201 lb . which were estimated using the 
half-mass radius and global velocit y dispersion ctoo- The 
o"p,obs values from lStrader et al.l (1201 1 ) were the central veloc- 
ity dispersions estimated using the observed velocity disper- 
sions by integrating a known King model over the fiber aper- 
ture, and are consistent with the predicte d line-of-sight veloc- 
ity di spersions at the cluster center from IBarmbv et all (120071 
120091) and this paper However, there are few young clusters 
with velocity dispersion measured by Strader et al. (201 1 ), so 
more observation data and analysis are needed to check the 
conclusion of the comparisons. 

Figure 6 plots the correlations of velocity dispersion and 
mass with M/L for M3 1 clusters. The left two panel s show 
these parameters of clusters from IStrader et al.l (t201 lb . while 
the right two panels show those derived by King model from 
B armby et al. (2007, 2009) and this paper We can see that 
the CTpfl and Mmod from King-model fits show large depen- 
dence on the MIL. The old and young clusters show a dis- 
tinct boundary of the M/L va lues. There is no cle ar correla- 
tion for MIL and cTp ^hs from' Strader et al.l (1201 lb . while the 
decrease in M/L toward lower masses is expected from dy- 
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Fig. 6. — Correlations of velocity d ispersion and mass with M/L in V 
band. Left panels show clusters from Strader et alj J2011I) (a sterisks) while 
righ t panels show clusters from Barmbv et al. (2007^ (crossesl lBarmbv et al.l 
12009ft (triangles), and this paper (filled circles). 

namical evolution like evaporation, which is the steady loss 
of low-mas s (high M/L) stars from the cluster driven by re 
laxation (Porteg ies et al.ll2bl0i) 



pendence on M/ L of cr^ 



Considering the sensitive de- 
■pfl, in the following discussion of the 
correlations of velocity dispersion with ellipt icity and [ Fe/H1 , 
we decided to use the CTn nhs values from :Strader et al.l (1201 lb 
for clusters in lBarmby et al.l (120071 120091) and this paper. 

3.4. Comparison of Three Model Fittings 

In order to deter mine which model can describe the struc- 
ture of clusters best. iMclaughlin & van der Marell (12005 ) and 
iMclaughUn etall (120081) defined a relative index, which 
compares the of the best fit of an "alternate" model with 
the x^ of the best fit of King model, 

A = OTaltemate ' xiei) I ^l\i^rn:i, +xi(,6)- (10) 

It is evident that the "alternate" model is a better fit than 
King model if A is negative, while King model is a better fit 
if A is positive. 

Figure 7 shows the relative quality of fit, A for Wilson- and 
Sersic-model fits versus King-model fits for the sample clus- 
ters in this paper The A values are plotted as a function of age 
and some structures, including the half-light radius 7?/„ the to- 
tal model luminosity Lmod, and the surface brightness over the 
half-light radius in the V band < fiy >/,. The circles refer to 
clusters with /?iast > 5/?/,, while triangles present those with 
/?iast < 5/?/,, where /?iast is the most large radius for the avail- 
able observation data. It can be seen that most clusters with 
/?iast < 5/?/, are those ha ying fainter luminosity Lmnd or sur- 
face brightness < jUy >h- IMclaughlin & van der Marell (l2005l) 
presented that the fitting data out to T^iast > 5Rh can effectively 
constrain the model fitting to the outer regions, which are es- 
sential to determine which structural model best describes the 
clusters. However, when /?iast is small (< 5Rh), few data are 
available at large cluster radii, and the fitting are mostly de- 
pendent on the inner part. So, these clusters cannot be used 
to determine a preference of one model or the other Simi- 
larly, we cannot conclude that these clusters are well fitted by 
both two models even A for them are small. .Mclaughlin et al.. 
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Fig. 7. — Relative quality of fit for Wilson and Sersic models (grey and open 
circles refer to the clusters with Ri^st > 5R;,, while grey and open triangles 
the clusters with Sjast < 5i?;,) versus King models for the sample clusters in 
this paper. 



(I2OO8I) determined a catalogue of structural and dynamical pa- 
rameters for GCs in NGC 5128 using the three models, and 
showed that the bright clusters (Lmod > 10^ Lq) are better fit- 
ted by Wilson model and Sersic model, indicating that the 
halos of clusters in NGC 5128 are more extended than what 
King model describes. There is no evident correlation of A 
with L mod and < >/, for clusters in this paper. Several 
studies dElson. Fall & Freemanll987l;ISan Rornan et al.l2012l) 
showed that young clusters in the Large Magellanic Cloud 
(LMC) and M33 do not appear to be tidally truncated and 
seem to be better fitted by power-law profiles than King mod- 
els, while old clusters show no clear differences between the 
qualities of the fittings. However, there is no correlation of A 
with age in Figure 7. We should notice that there are only six 
young clusters in our sample, and 13 clusters with no ages es- 
timated in previous studies are assumed to be 13 Gyr A large 
sample of young star clusters with precise age estimates are 
needed for the study on correlation of A with age. However, 
we do see that the King- and Wilson-model fits are better than 
the Sersic-model fits. We concluded that clusters in M31 can 
be well fitted by both King model and Wilson model (also 
reported in Barmby et al. 2007, 2009). 

Figure 8 compares the relative quality of fit, A values with 
a number of structure parameters [Re, R/,, l^vfi, i^mod, cTp.o, 
and Eb) for the sample clusters in this paper The grey and 
open circles show the physical properties of clusters with 
^last ^ 5/?/, derived from Wilson- and Sersic-model fits com- 
paring to King-model fits, respectively. The triangles refer 
to clusters with /?iast < 5/?/,. There are some clusters with 
comparable x^^ but large discrepancy of Rf, and L^od val- 
ues for King- and Wilson-model fits. We can see that most 
of these clusters have /?iast < 5/?/,. As discussed above, 
the few constrain of the fitting to outer regions results in 
much different extrapolations of models, and it is hard to de- 
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Fig. 8. — Comparison of some parameters for Wilson and Sersic models 
versus King models for the sample clusters in this paper, including the pro- 
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brightness in the V band fiv,o, the total model luminosity Ln,odi the predicted 
central line-of-sight velocity dispersion CTpo, and the global binding energy 
Ej,. Symbols are as in Fig. 7. 

termine which model does th e correct fitting in those cases 
dMclaughhn & van der Marell [2005). Most parameters de- 
rived from Wilson model are slightly larger than those from 
King model, while parameters derived by Sersic model are 
smaller than those from Kin g model, especially for R^, cTpfi, 
and Eh- iBarmbv et al.l (120071) concluded that rather than an in- 
trinsic difference between clusters in M3 1 and other galaxies, 
the preference for King models over Wilson models for M3 1 
clusters is due to some subtle features of the observations. 

4. DISCUSSION 

We combined the newly derived parameters here with those 
derived b y King-mo del fits for M3 1 young massive clusters 
(YMC s) fearmby et al. 2009), M31 globulars (Barmbv et al. 
120071) . M31 extende d clusters (ECs) (.Huxor et al.. ,201 11 
and MW globulars dMclaughhn & van der Marell 120051) to 
construct a large sample to look into the correlations be- 
tween the parameters. The ellipticities and galactocentric 
distances for the MW GCs are from Harris) (119961) (2010 
edition). The p arameters used in th e following discussion 
for M31 GCs (iBarmbv et all l2007h are those derived on 
the bandpass close to V band (e.g., F555W and F606W), 
w hile the data fr o m blu er filters are preferred for the YMCs 
in iBarm by et al.' (2009"). The metallicities for most young 
clusters in Barniby et al. (2009) were assumed to be solar 
metallicity, while metallicities from Perina et al. (2010) were 
adopted for s ome older clusters (B083, B222, B347, B374, 
and NB16). iHuxor et"an d201 lb derived structure parame- 
ters of 13 ECs, including the core radius, half-light radius, 
tidal radius, and the central surface brightness in magnitude. 



The metallicities of 4 ECs (HEC4, HEC5 HEC7 and HEC12) 
were determined by iMackev et al.l (120061) . and the integrated 
cluster masses for them were derived here using the abso- 
lute magnitudes and M/L values in V ba nd from population- 
synthesis models (see lWang & Mall2012L for the details). 

4.1. Galactocentric Distance 

Figure 9 shows structural parameters as a function of galac- 
tocentric distance R^^ for the new large sample clusters in 
M31 and the MW. Some global trends can be seen. Both R/, 
and r, increase with the R^^ as e xpected, although t he trend is 
largely driven by the ECs from lHuxor et al.l (1201 lb . The star 
clusters with large Rg^ can keep more stars with weake r tides 
relatively to those located at bulge and disk. Georgiev eLalJ 
(2009) explained that the Rf, and r, of the clusters in the halo 
of galaxies, which may be accreted into the galaxies from 
dwarf galaxies, can expand due to the ch ange from strong to 
a weaker tidal field. Barmbv e t al.l (120021) concluded that the 
correlation of R/, with R„c refle cts physical f ormation condi- 
tions as suggested by van den Berghl d 19911) for GCs in the 
MW. However, S trader et al. (20,1 2) found that no clear cor- 
relation between R/, and Rgc exists beyond Rg^ ~ 15 kpc for 
GCs in NGC 4649, and they suggested that the sizes of GCs 
are not generically set by tidal limitation. It can be seen that 
there are no compact clusters at large radii (/?gc > 40 kpc) in 
the MW, while in M31, there is a bimodality in the size dis- 
tribut i on of GCs at large radii (also reported in iHuxor et"aLl 
1201 U IWang & M"all2012h . It is interesting that there are few 
GCs having R^ in the range from 8 to 15 pc at large radii 
(Rgc > 40 kpc) in M31. There are three clusters (B124, 
B127, and NB39) with large R/, and r, at small Rg^, leading 
to more diffuse trends. However, the total model masses for 
them are 10^'^^, 10^'^^, and IQ^ '^^Mq, respectively, meaning 
that they are massive clusters. It is not unexpected that these 
massive clusters can contain more stars than other clusters, al- 
though they are at small Rg^. The luminosities of clusters de- 
crease with increasing Rg^, implying that either strong dynam- 
ical friction drives predominantly more massive GCs inwards, 
or massive GCs may favor to form in the nu clear regions of 
galax ies with the higher pressure and density (iGeorgiev et alJ 
l2009h . However, the lack of faint clusters with small galac- 
tocentric distance may be due to selection effects, since these 
clusters are difficult to detect against the bright back ground 
near M 31 center, which is also reported by IBarmbv et al.l 
d2007h using the correlation of central surface brightness with 
galactocentric distance. The metallicities of clusters decrease 
with increasing Rg^, indicating that metal-rich clusters are typ- 
ically located at smaller galactocentric radii than metal-poor 
ones, although with large scatter. Ivan den Berghl (11991) found 
that metal-rich clusters ([Fe/H]> -0.8) in M31 seem to form 
in a rotating disk extending to Rg^ ^ 5 kpc. The metal-rich 
GCs may have underg one accelerated inter nal evolution due 
to strong tidal shocks dStrader et al.l 1201 ih from both bulge 
and disk passages. There are two matal-poor clusters (Bl 14, 
B324) at small Rg^, which are at the bottom-left in the [Fe/H]- 
Rgc panel. The diffuse distribution of these parameters vs. Rg^ 
may be caused by the different galactocentric distance used, 
which is true three-dimensional distances for Galactic GCs 
while projected radii for M31 clusters. 

4.2. Ellipticity 

Figure 10 shows the distribution of ellipticity with galacto- 
centric distance, metalhcity, and some structure parameters 
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the observed velocity dispersion, and stracture parameters. Symbols are as in 
Fig. 9. 

for clusters in the MW and M31, which may show clues 
to the primary factor for the elongation of clusters: rota- 
tion and velocity anisotropy, cluster mergers, "remnant elon- 
gation", and galaxy tides (Larsenetal. 2001; Barmby et al. 
120071) . iGeyer et alJ (l2009i) presented that the tidal forces can 
only distort a cluster's outer region, while the elongation of 
the inner parts are due to its own gravitational potential and 
the total orbital angular momentum. 

1) Rotation and velocity anisotropy. Cluster rotation 
is the g enerally accepted explanation for cluster flattening 
(iDavous t & Pmgniel 1990). Dynamical models show that in- 
ternal relaxation coupled to the external tides will drive a 
clus ter toward a round er shape over several re laxation times 

i see [ Harris et al.ll2002l and references therein). iBarmbv et alj 
2002!) presented that more compact clusters which expe- 
rience more relaxation processes, and clusters with larger 



velocity dispersions which r otate more slowly, should be 
rounder. iBarmbv et al.l (12007b presented that dynamical evo- 
lution could reduce the initial flattening caused by rotation 
or velocity anisotropy, indicating that more evolved clus- 
ters would be rounder In order to check these predic- 
tions, we show the correlations of ro, c, the observed ve- 
locity dispersion cr^^obs, and t^ji with the ellipticity. The ob- 
served veloci ty disp ersion crp_obs for M31 clusters are from 
[Strader et al l (1201 lb. while cTp pt^^ for clu sters in the MW are 
from Mclaughli n & van der Marell (120051) . All the o-^ obs have 
been extrapolated to their central values with an aperture of 0. 
There is no clear correlation of ellipticity with these parame- 
ters. However, we do see that more compact clusters-smaller 
ro and larger c- are more rounder, which is consistent with 
previous studies dBarmbv et al.ll2002i) . 

2) Cluster mergers and "remnant elongation". 
Ivan den Bergh & Morbeyl ([1983) presented that elliptic- 
ity correlates strongly with luminosity for clusters in 
LMC: more luminou s clusters are more flattened, while 
Ivan den Berghl (119961) presented that the most flattened GCs 
in both the MW and M31 are also brightest. This may be 
explained by the cluster mergers and "remnant elongation", 
which is from some clusters' former lives as dwarf galaxy 
nuclei. So, the larger ellipticities for clusters located at large 
galactocentric radii of M3 1 may be related to t he merger or 
accretion history tha t M3 1 has experien ced (e.g. lBekkill20Tol : 
iMackev et all 1201 Ot iHuxor et all I2OI V). However, no clear 
correlation of ellipticity with luminosity can be seen. 

3) Galaxy tides. A strong tidal field might rapidly de- 
stroy the velocity anisotropics, and force an initially triax- 
ial, rapidly rotating elliptical GC to a more isotropic distri- 
bution and spherical shape, while weak tidal fields are un- 
able to change the initial shapes of GCs (lGoodwinlll997h . It 
seems plausible that clusters located at diff'erent galactocen- 
tric radii or different galaxies have variou s distribution s of el- 
lipticities, due to the dive rse galaxy tides. iHarris et"aD (120021) 
and Ba rmbv et al.l (l2007h found the distributions of elliptici- 
ties for M31 and NGC 5128 are very similar, but differ from 
the MW distribution, which has few very round clusters. No 
clear correlation of elliptici ty with R^c. can be s een for clusters 
in these galaxies (also see Barmby et al."'2007). Howe ver, the 
innermost clusters are slightly more spherical (B armby et alJ 
l2002h . which may due to the strong tidal field near galaxy 
center that reduces ellipticities of them. Some clusters lo- 
cated at large projected radii of M31 do have larger elliptici- 
ties than most GCs in M3 1 and the MW, which may be caused 
by galaxy tides coming from satellite galaxies (Wang & M^ 
2012). 

San Roman et alj ( 120121) summarized the orientation of a 
large sample of clusters in M33 and found that the distribu- 
tion of PAs shows a strong peak at -55°, which is close to 
the direction towards M3 1 . These authors suggested that the 
elongation of clusters in M33 may be attributed to the tidal 
forces of M31, considering that a recent encounter between 
M33 andM31 ([M cConnachie et al. 2009i jPutman et alj|2009l: 
ISan Roman et al.ii2010i:iBernard et al.n2012l) mav have led to 
significant effects on the properties of M33 disk. Figure 1 1 
depicts the distribution of PAs of star clusters in this paper 
and 34 clusters from Barmbv et al. (2007), which were not 
included in Barmby et al. (2002). No clear trend is present in 
the orientation vectors towards M33, although a small peak 
at -60° does exist in the distribution of the PAs. Here we 
concluded that the elongation of clusters seems to be due to 
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Fig. 1 1 . — (Left panel) Cluster elongations and orientations shown with 
respect to M31. The vector sizes are correlated with the ellipticities. The 
ellipse has a radius of 10 kpc and is flattened to b/a = 0.6. (Right panel) 
Distribution of the P.A.s of star clusters in M31. 

various factors (iBarmbv et al.ll200'7l) . 

4.3. Metallicity 

Figure 12 plots structural parameters as a function of [Fe/H] 
for the new large sample clusters in M3 1 and the MW. The 
trends of the parameters with [Fe/H] nearly disappear when 
we add the data f or the M3 1 YMCs, which were fitted with 
solar metallicities ('Barm bv et al.ll2009h . However, we noticed 
that the me tallicities for most of these YMCs obtained by 
iKang et alJ(l2012h are poorer th an [Fe/H]= -1.1. If we do not 
consider the YMCs in lBarmby et al.. (.2009) . it is evident that 
the metal-rich clusters have smaller average va l ues of Rh than 
those of metal-poor ones (iLarsen et al.l 120011: iBarmbv et al.l 
120021) . IStrader et af] (120 12h presented that the sizes of metal- 
rich GCs are smaller than the metal-poor ones in NGC 4649, 
which is a massive elliptical galaxy in the Virgo galaxy clus- 
ter These authors explained th at as an intrinsic siz e difference 
rather than projection effects. Sippel etaP (120121) carried out 
N-body simulations of metallicity effects on cluster evolution, 
and found that there is no evident difference for the half-mass 
radii of metal-rich and metal-poor cluster models. So, they 
explained that metal-rich and metal-poor clusters have similar 
structures, while metallicity effects combined with dynamical 
effects such as mass segregation produce a larger difference 
of the half-light radii. 'Barmby etaP (12007*) also found that Ri, 
decreases with increased metalhcity for GCs in the MW, the 
MCs, Fornax dwarf spheroidal, and M31, except for GCs in 
NGC 5128. An evident feature is that these four ECs are all 
metal-poor and have large /?/,. No clear correlation of model 
luminosity Lmod with [Fe/H] is present, but we do see that the 
metal-rich clusters tend to be more luminous. Both the ob- 
served velocity dispersion o-p^obs and central "escape" veloc- 
ity Vesc.o increase with metallicity. Georgiev et al. (2009) pre- 
sented that the higher Vesc of more metal-rich clusters in the 
observation may reflect the metallicity dependence of the ter- 
minal velocities of the stellar winds, since the Vesc of a metal- 
rich cluster should be higher to retain the fast stellar winds. 

4.4. Core-collapsed Clusters 

Core-collapsed clusters in general show a power-law slope 
in the central surface brightness profiles, which can be better 
fitted by a power-law model than King model (Barmb y" et alJ 
|2002). Noyola & Gebhai'dt (2006) presented that the process 
of core collapse can be divided into two stages. In the first 
stage, stars are driven to the halo of the cluster due to close 
encounters. Stellar evaporation occurs and the core contracts 
due to energy conservation. In the second stage, the low- 
mass stars are scattered to high velocities and escape to the 
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Fig. 12. — Structural parameters and the observed velocity dispersion as a 
function of [Fe/H]. Symbols are as in Fig. 9. 

halo, while the high-mass stars sink to the core due to energy 
equip artition. The increasing core density also increases the 
interaction rate of the binaries and single stars, which can gen- 
erate energy in the core, reverse the contraction process, and 
produce an expansion. After a long time, the core shrinks 
again, and the process repeats. M15, which has a central 
cusp, may be at an interm ediate state betw e en the extremes 
of coll apse and expansion (iDuU et alj|1997l) . IChatteriee et al.l 
(120121) presented that the core-collapsed clusters are those that 
have reached or are about to reach the "binary burning" stage, 
while the non core-collapsed clusters are still contracting un- 
der two-body relaxation. Trager et al. (1995) presented a cat- 
alogue of surface brightness profiles of 125 Galactic GCs, and 
classified 16% of their sample a s core-collapsed clusters and 
6% as core-collapsed candidates. iMclaughlin et al] ( 120081) no- 
ticed that a number of clusters (more than half of their sample) 
in NGC 5 128 have the index n > 2 derived from Sersic model, 
indicating that these clusters have strongly peaked central 
density profiles. These authors explained that the PSF may 
flatten the models, and then the cuspy centr a l dens ity profiles 
are shown relatively. iNovola & GebhardB (12006b presented 
that Galactic clusters with steep cusps are all close to the cen- 
ter of the Galaxy, indicating that an increased incidence of 
tidal shocking (Gnedin et aL 19 99) may accelerate the core 
collapse process. Barm by et al.l (^002) also concluded that 
most of the M3 1 core-collapsed candidates are within 2 kpc 
of the center of M3 1 . 

Figure 13 shows the distributions of galactocentric distance, 
and some structure parameters derived from Sersic-model fit- 
ting for clusters in this paper and M31 GCs in Barmby et al] 
(12007) . There are 25 clusters bested fitted by Sersic model 
^ ^K66 ^"'^ ^ ^w^' ™^ 1 1 of them have n > 2, 
indicating that they have cuspy central density profiles. We 
assumed these 1 1 clusters to be possible core-collapsed clus- 
ters. In fact, the clusters with n > 2 do have smaller Rr and 
larger /?/,/^c than their counterparts with n < 2. Barmbvet al.l 
(2002) presented several clusters (BOll, B064, B092,"B123, 
B145, B231, B268, B343, and BH18) as core-collapsed can- 
didates in M3 1 . However, most of these clusters are best fitted 
by King model or Wilson model in this paper except B 145 and 
B343. We checked the images of these clusters and found that 
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1 10 100 1 10 100 1000 10^ 10^ 



Fig. 13. — The distributions of the galactocentric distance and several stiuc- 
ture parameters derived from Sersic- model fitting for clusters in this paper 
and M3 1 GCs in lBarmbv et"J] ROpTh . 

the images of B123, B231, and B268 are not fully resolved, 
while the SNRs of images for B064 and B092 are low. Ac- 
tually, the discrimination between the core-collapsed clusters 
and the "King model clusters" (Mevlan & Heggie 1997) may 
often become unclear for several reasons: (1) statistical noise 
due to some unresolved bright stars in the cores of clusters; 
(2) the similarity between th e high-concentration King -model 
and power-law profiles (see iMevlan & Heggielll997l for de- 
tails). B145 is best fitted by Sersic model in this paper, but 
n - 1.5. The cluster B343 shows a cuspy density profile in 
the center, and has been classified as a core-collapsed cluster 
dBendinelH et al. 1993; Grillmairet al. 1996). We also find a 
better fitting by Sersic model for it with n - 3.35. The distri- 
bution of core-collapsed cluster candidates is not constrained 
to the center of M3 1 . There are only two out of the 1 1 clusters 
with /?gc < 2 kpc, while four ones with R^^ > 15 kpc. 

4.5. The Fundamental Plane 

It has been widely noticed that GCs do not occupy the full 
four-dimensional parameter space (concentration c, scale ra- 
dius ro, central surface brightness fiyfl, and central M/L or 
velocity dispersion ctq) instead locate in a remarkably nar- 
row "fundamental plane" (FP). It is interesting to learn which 
of the structural and dynamical properties are either univer- 
sal or dependent on characteristics of their parent galaxies 
(^eylan & Heggie 1997). Djorgovski & Meylan ( 1994) ex- 
plored some tight correlations between various propert ies us- 
ing a large sample of Galactic clusters, and Djorgovskii (119951) 
defined an FP for GCs in ter ms of velocit y dispersion, ra- 
dius, and surface brightness. ISaitol (119791) investigated the 
mass-bindi ng energy relations hip for a few GCs and ellipti- 
cals, while McLaughlinI (120001) defined an FP using the bind- 
ing energy and luminosity , which is fo rmally different from 
that of Djorgovski ( 1995). iHarris et al.l (2002) found that the 
NGC 5128 GCs describe a relation between binding energy 
and luminosity tighter than in the MW. Here we show the two 
forms of the FP for the new large sample clusters in the MW 
and M31. 




Eq/Mq pc 2 s^/Mg pc 




10^ 10^ 10^ 10"* 10^ 10^ 10^ 10-^ 10* 10^ 



<S>/Mq pc 2 <^>h/^@ ^ 

Fig. 14. — Evidence of an FP of the cluster parameters, which is defined in 
terms of central velocity dispersion CTp o, radius or R/,, and surface mass 
density < Z >o or < Z >;,. Symbols are as in Fig. 9. 

Figure 14 shows the mass-density-based FP relations. The 
left two panels show the correlations of the properties de- 
rived with the M/ L values from the velocity dispersions 
(lStraderetal.ll20m . using King model and the same fitting 
process in this paper (Section |3]l, while the right two pan- 
els show the correlations of the properties de rived with the 
M/L v alues from population-synthesis models. iBarmbv et al.l 
( 120091) presented the surface-brightness-based FP relations 
and found a large offset between the young M3 1 clusters and 
old clusters. They explained that as a result from lower M/L 
values for the young clusters in M3 1 . It is obvious that the 
velocity dispersion, characteristic radii, and surface density 
for these clusters show tight relations, both on the core and 
half-light scales. The exist of FP for clusters strongly reflects 
some universal physical conditions and processes of cluster 
formation. 

Figure 15 shows the correlation of binding energy with the 
total model mass. The left panel shows the Mmod and Ei, de- 
rived wit h the M/L from th e directly observed velocity dis- 
persions (i StraderetaLllMTl) . using King model and the same 
fitting process in this paper (Section |3}, while the right panel 
shows those properties derived with the M/L from population- 
synthesis models. All the clusters locate in a remarkably tight 
region although in the widely different galaxy environ ments, 
which is co ns istent with the previous studies (Barmb v et al.l 
120071 120091) ■ IBarmbv et al] (|2007) concluded that the scat- 
ter around this relation is so small that the structures of star 
clusters may be far simpler than those scenarios derived from 
theoretical arguments. 

5. SUMMARY 

High-resolution imaging can be derived from //ST observa- 
tions for M3 1 star clusters. In this paper, we presented surface 
brightness profiles for 79 clusters, which were selected from 
.Barmbv et al..(.2002) and .Mackev et al..(.2007) . Structural and 
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10^ 10^ 10^ 10^ 10* 10^ 10*' 10^ 



mod' mod' 

Fig. 15. — Evidence of an FP of the cluster parameters, which is defined in 
terms of binding energy £4, and the total model mass Mmod ■ Symbols are as 
in Fig. 9. 



and the Galactic clusters fall in the same regions of param- 
eter spaces. There is a bimodality in the size distribution 
of M31 clusters at large radii, which is different from their 
Galactic counterparts. There are 1 1 clusters in M31 best fitted 
by Sersic models with index n > 2, meaning that they have 
cuspy central density profiles, which are classified as core- 
collapsed cluster candidates. We investigated two forms of 
the FP, including the correlation of velocity dispersion, radius, 
and surface density, and the correlation of binding energy 
with the total model mass. The tight correlations of cluster 
properties indicate a tight FP for clusters, regardless of their 
hos t environments, which is c onsistent with previous stud- 
ies jBarmbv et al.ll2007l 120091) . In addition, the tightness of 
the relations for the internal properties indicates some physi- 
cal conditions and processes of cluster formation in different 
galaxies. 



dynamical parameters were derived by fitting the profiles to 
three different models, including King model, Wilson model, 
and Sersic model. We found that in the majority of cases. 
King models fit the M31 clusters as well as Wilson models, 
and better than Sersic models. 

We discussed the properties of the sample GCs her e com- 
bined with GCs in the MW (Mclaughlin & van der Mare| 
20051) and clusters in M31 (Barmbvetal. 2007, 2009; 
Huxor et all |20TT1) . In general, the properties of the M31 
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TABLE 1 

Integrated measurements for 79 sample clusters em M31. 



Name 




ep 




0^ 


B 


V 


7 




E{B - V) 


[Fe/H] 


Age 






(degEofN) 




(degEofN) 


(Vegamag) 


(Vegamag) 


(Vegamag) 


(kpc) 






(Gyr) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


B006 


0.14 + 0.01 


-40 ±8 


0.24 ± 0.02 


-45 ±8 


16.49 


15.50 


14.33 


6.43 


0.11 


-0.59 ± 0.41 


13.25 


BOll 


0.12 ±0.01 


21 ±8 


0.15 ±0.01 


32 ±7 


17.46 


16.58 


15.56 


7.73 


0.09 


-1.71 ±0.24 


15.85 


B012 


0.17 ±0.02 


24 ±8 


0.14 + 0.01 


18 ±7 


15.86 


15.09 


14.03 


5.78 


0.11 


-1.91 ±0.21 


8.00 


B018 


0.19 ±0.01 


-9 ±5 


0.15 ±0.02 


-5 ±8 


18.37 


17.53 


16.38 


9.31 


0.20 


-0.77 ± 0.39 


1.20 


B027 


0.18 ±0.02 


11 ±7 


0.16 ±0.03 


-12 ±9 


16.49 


15.56 


14.41 


6.02 


0.19 


-1.64 ±0.16 


15.75 


B030 


0.16 ±0.01 


-12± 12 


0.13 ±0.01 


-37 ±6 


18.75 


17.39 


15.68 


5.66 


0.57 


-0.14 ±0.26 


9.62 


B045 


0.12 ±0.01 


-13±6 


0.16 ±0.02 


-13±7 


16.72 


15.78 


14.51 


4.90 


0.18 


-1.01 ±0.50 


11.40 


B058 


0.09 + 0.01 


-27 + 8 


0.11+0.01 


-36 + 9 


15.81 


14.97 


13.87 


6.96 


0.12 


-1.02 ±0.21 


8.01 


B068 


0.21 + 0.02 


31+3 


0.18 ±0.02 


37 + 4 


17.60 


16.39 


14.84 


4.32 


0.42 


-0.41 ±0.17 


7.90 


B070 


0.12 + 0.01 


1 ± 10 


0.16 + 0.02 


-24 ±11 


17.61 


16.76 


15.72 


2.46 


0.12 


-1.42 ±0.43 


8.75 



and 6 of bluer filters, 
and d of redder filters. 
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TABLE 2 
CAUBRAnoN Data for HST images. 



Filter 


Pivot A 
(A) 


«5 


Zeropoint* 




Conversion Factor^ 


Coefficient" 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


Calibration Data for ACS/WFC images 


F435W 


4318.9 


4.20 


25.779 


5.459 


3.1693 


27.031 


F475W 


4746.9 


3.72 


26.168 


5.167 


1.6926 


26.739 


F555W 


5361.0 


3.19 


25.724 


4.820 


1.8508 


26.392 


F606W 


5921.1 


2.85 


26.398 


4.611 


0.8207 


26.183 


F814W 


8057.0 


1.83 


25.501 


4.066 


1.1349 


25.638 


Calibration Data for WFPC2 images 


F300W 


2986.8 


5.66 


21.448 


6.061 


297.9614 


27.633 


F450W 


4555.4 


3.93 


24.046 


5.263 


13.0545 


26.835 


F555W 


5439.0 


3.14 


24.596 


4.804 


5.1542 


26.376 


F606W 


5996.8 


2.81 


24.957 


4.581 


3.0100 


26.153 


F814W 


8012.2 


1.85 


23.677 


4.074 


6.1342 


25.646 



=RaX E(B - V). 

''Additive conversion between surface brightness in counts s"' arcsec"^ and magnitude in mag arcsec"^. 
■^Updated absolute magnitude of the sun (C. Willmer, private communication). 

''Multiplicative conversion between surface brightness in counts s"' arcsec"^ and intensity in L© pc"^. 
■^Additive conversion between surface brightness in magnitude in mag arcsec"^ and intensity in pc"^ . 



TABLE 3 

The iNTENSrrY profiles for 79 sample clusters in M31. 



Name 


Detector 


Filter 


R 


/ 


Uncertainty 


Flag 








(arcsec) 


(Lo pc-2) 


(Lo pc-2) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(V) 


B006 


WFPC2/PC 


F555W 


0.0234 


44001.773 


155.603 


OK 




WFPC2/PC 


F555W 


0.0258 


43776.562 


171.971 


DEP 




WFPC2/PC 


F555W 


0.0284 


43528.629 


188.743 


DEP 




WFPC2/PC 


F555W 


0.0312 


43258.270 


206.935 


DEP 




WFPC2/PC 


F555W 


0.0343 


42964.875 


227.137 


DEP 




WFPC2/PC 


F555W 


0.0378 


42647.070 


249.659 


DEP 




WFPC2/PC 


F555W 


0.0415 


42303.609 


274.847 


DEP 




WFPC2/PC 


F555W 


0.0457 


41919.195 


299.011 


DEP 




WFPC2/PC 


F555W 


0.0503 


41381.914 


308.833 


OK 




WFPC2/PC 


F555W 


0.0553 


40666.176 


309.009 


DEP 




WFPC2/PC 


F555W 


0.0608 


39884.902 


322.641 


DEP 




WFPC2/PC 


F555W 


0.0669 


38966.234 


328.825 


DEP 




WFPC2/PC 


F555W 


0.0736 


37925.949 


316.808 


OK 



TABLE 4 
Coefficients for the PSF models. 



Detector 


Filter 


ro 

(arcsec) 


a 


j8 


(1) 


(2) 


(3) 


(4) 


(5) 


ACSAVFC 


F435W 


0.068 


3 


3.80 




F475W 


0.064 


3 


3.60 




F555W 


0.057 


3 


3.39 




F606W 


0.053 


3 


3.14 




F814W 


0.056 


3 


3.05 


WFPC2/WFC 


F300W 


0.076 


2 


5.05 




F450W 


0.073 


2 


4.89 




F555W 


0.064 


2 


4.35 




F606W 


0.059 


2 


4.11 




F814W 


0.051 


2 


3.71 


WFPC2/PC 


F300W 


0.051 


2 


3.76 




F555W 


0.045 


2 


3.10 




F606W 


0.045 


2 


2.96 




F814W 


0.059 


2 


3.18 



TABLES 

Basic parameters of 79 sample clusters m M31. 



Name Detector Band Nj.t," Model atL* ^bkg' Wo'' cjtf ix^' logro' logro* 

(Lo pc"^) (mag arcsec"^) (arcsec) (pc) 
(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) (12) 



B006 


WFPC2/PC 


555 


49 


K66 


35.44 


5.10+1.54 


7 90+°-"* 


1.80!«;JJ^ 


1 4 4^+0.04 


-0.900+°:°}^ 


-0 320+°'"- 








49 


W 


5.53 


-8.30 ± 0.92 


7-60!°:°^ 


2.86!S!;S^ 


14 44+0-01 


-0.850!°:°» 


-0 270+" °°' 

"•^'"-0.007 








49 


S 


130.08 


3.50 ±2.56 


8.20!«:^^ 


2 60+° °* 

^•""-0.04 


14.28!«:i^ 


-2 400+'' ''^2 




B006 


WFPC2/PC 


814 


49 


K66 


14.97 


6.00 ± 3.36 


1.89!°- 


13 35+°-''* 


-0.950!«:°» 


-0 370+°°'^ 

"•■'"-'-0.026 








49 


W 


1.88 


-9.90 ± 2.25 




2.861S^ 


1 ^ QQ+0.04 


-o.850!S:°S? 


-0 270+° °°" 

"•^' -0.007 








49 


S 


32.24 


4.40 ± 3.68 






16 44+0-50 


-2 4sn+0-038 


1 8^n-^0-038 

-i-»/0_nmq 



"The number of points in tlie intensity profile that were used for constraining the model fits. 
''The minimum obtained in the fits. 
■^The best-fitted background intensity. 

''The dimensionless central potential of the best-fitting model, defined as Wo = -^(Q)/cr^. 
"The concentration c = log(r,/ro). 

^The best-fit central surface brightness in the native bandpass of the data. 
SThe best model-fit scale radius tq in arcseconds. 
'"The best model-fit scale radius ro in parsecs. 



On 



TABLE 6 

Derived structural and photometric parameters of 79 sample clusters in M31. 

Name Detector Band Model lograd" log^T iog^? log/? J/?/ logk' logTT logLy« log 4' < Mv >h^ 

(pc) (pc) (pc) (Ley pc-2) (Lo.v pc"') (Lq.v) (mag) (Lg.,/ pc'^) (mag arcsec"^) 

(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) (12) (13) (14) 

B006 WFPC2/PC 555 K66 1.48+«-«5 -0.333!«-«|« 0.386!»!M« 0.719!««^^ 4.79!««3 4.81!«-«5 5.63!««2 l5.18+«-«^ 4.06t_°i^ 16.20!«-23 

W 950+0.04 _n9q<:+0.008 n<:il+0.046 n Rn7+''053 4 70+O.O2 4 77+O.O3 <:71+0.02 UQQ+O.OS T Sq+014 if,f,o+0.17 

yv u.z,7j_QpQ.^ "•-'^'^-0.086 "-""'-O.OgS -0.02 ^-''-0.03 -'•"^-O.Oa •'^°-0.06 -0.07 ^"•"-'-0.35 

S cx, -1.820tg3? 0.399tJ^« 2.218^;?^ 4.85^^ 5.49tg^ 5.63t«^ 15.18t!^ 4.04lJ« 16.27t« 

B006 WFPC2/PC 814 K66 1.52+;;-o« -0.38r;;-0ii 0.428+o;;«^ 0-810+!Ji^^ 4.79+«;;3 4.86+;;-»« 5.63+o;;2 15.18+;;-»^ 3.98+»[^ 16.41+0-^J ^ 

W 7 *;q+o:o9 _o oqj: +0:005 q^i -1+0:093 n 8nQ+'':o98 470+0:02 477+0:03 £571+0:02 .400+0:08 o oq+o:o6 if:f;4+o:41 g 

V¥ /..js- Qpg u.z,7J_QpQj U.JlJ pQ^^ U.OU7_p(^g ^•'°-0.02 ^-"-0.03 -'•"-O.OS '^"-O.Oe ■'•°^-0.16 ^"-^ -0.14 3 

S cx, -1.870!»:°3B 0.416!°;}^^ 2.286t;^^ 4.85t^ 5.53!»:°^ 5.63^^ 15.18!»:i^ 4.00::°;}« 16.35!»:^^ ^ 

"The model tidal radius r, in parsecs. ^ 
''The projected core radius of the model fitting a cluster, defined as I(Rc) = /o/2. 

'^The projected half-light, or effective, radius of a model, containing half the total luminosity in projection. 
''a measure of cluster concentration and relatively more model-independent than Wo or c. 

'The best-fit central (R = 0) luminosity surface density in the V band, defined as log/o = 0.4(26.358 - yUy.o). where 26.358 is the "Coefiicient" corresponding to a solar absolute magnitude Mv,e = -1-4.786 (C. 
Willmer, private communication). 
'The central (r = 0) luminosity volume density in the V band. 
'^The V-band total integrated model luminosity. 

''The total V-hand magnitude of a model cluster, defined as Vtot = 4.786 - 2.5 log(Lv/Lo) + 5 log(Z)/10 pc). 

'The luminosity surface density averaged over the half-light/eft'ective radius in the V band, defined as log/^ = log(Lv/2nR^). 

JThe surface brightness in magnitude over the half-hght/elfective radius in the V band, defined as < jxy >h= 26.358 - 2.5 log 4. 



TABLE? 

Derived dynamical parameters of 79 sample clusters in M31. 



Name Detector Band T^" Model logM,o,* log^*" logSo'' logpo' logS,,^ logo-,;,o« logv^.^.o* log/,,,/ log/o^ 

(MoLgV) (^o) (erg) (Mepc-^) (Mq pc-^) (Mq pc-^) (kms-') (kms-') (yr) (Mo (pc km s^i)-') 5* 

(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) (12) (13) (14) >% 

B006 WFPC2/PC 555 2.588+^-«J^ K66 6.05+»|f 51.76+;;^« 5.20+^|J 523^, 4.48+;;-|^ ^■'^^2*"/>fl l-783+;;™'^ 9.08+^-;5 0.493+""*" ^ 

W 6.12-- 51.81-- 5.20!JJ:|? 5.19!«- 4.30^- 1.190!j;:«2 l-803!EiS 9.30^1^ 0.387!i|i S 

S 6.05!«;j^ 50.39-- 5.26^1? 5.90-- 4.45!^? 0.714!«:«« 1.480!":- 9.10!"- 3.054-;- ^ 

B006 WFPC2/PC 814 2.588-8^^ K66 6.05-;2 51.66-^^ 5.20+;!;2 5.28- 1] 4.39--2^ 1-148-|;S ^■'^^^Tofs 9-14-;^ OMT'^^lll 2 

W 6 12-->2 5181--'"* 5 20--'- 5 19— '" 430+0^20 1 loo+o"''" 1 803--'"''' 9 30--" 387--'"' 

S 6 05- " 50 Se- "*' 5 26+" '" 5 94- " 4 42+"-22 698+" °*° 1 470+" °'' 9 13+" "' 3 168+" °^' ^ 

"The V-band mass-to-light ratio. 

''The integrated cluster mass, estimated as log Mmt = log T^"'' + log Ly. ^ 
"^The integrated binding energy in ergs, defined as £j = -(1/2) ^' 4nr^p<f>dr. 

''The central surface mass density, estimated as log So = log T^"'' -1- log /q. S 

^The central volume density, estimated as log po = log T^°'' + log 70. £1 

^The surface mass density averaged over the half-light/effective radius Rh, estimated as log = log T^"'' -1- log 4. 
8The predicted line-of-sight velocity dispersion at the cluster center. 

•"The predicted central "escape" velocity with which a star can move out from the center of a cluster, defined as v'^scq/'^o ~ ^^^o + GMtot/r,ag]. 

'The two-body relaxation time at the model-projected half-mass radius, estimated as tr^ = i^Q^MlM/mt) '"In,'' ■ average stellar mass in a cluster, assumed to be 0.5Mq. 



JThe model's central phase-space density, defined as log/o = log|/Jo/(2?rcr^)""]. 



